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INTRODUCTION 
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Inventory  Smoke  Munitions  Test  (Smoke  Week  I)  was  held  at  Dugway  Proving 
Grounds  during  November  1977.  The  objectives  of  the  test  as  reported  in 
reference  a  were: 

1)  to  obtain  smoke  characteristics  of  inventory  smoke  munitions, 

2)  to  obtain  data  required  by  JTCG/ME  smoke  obscuration  model  for 
evaluation  of  existing  inventory  munitions, 

3)  to  obtain  attenuation  characteristics  of  inventory  smoke  munitions, 
and 

4)  to  obtain  data  on  the  effects  of  screening  smokes  on  the  optical 
properties  of  existing  and  developmental  electro-optical  systems. 

The  smoke  types  selected  for  analysis  from  Smoke  Week  I  are  HC  (hexachloroe- 
thane,  zinc  oxide,  aluminum  mixture),  WP  (white  phosphorus),  and  FWP  (foreign 
white  phosphorus).  Two  trials  each  of  HC  and  WP  were  conducted,  while  each 
FWP  was  tested  once.  A  summary  of  the  test  trials  appears  in  Table  1. 

An  aerosol  photometer  was  used  to  measure  the  concentration  loading  (CL)  or 
column  density  of  the  smoke,  while  various  transmissometers  measured  transmis¬ 
sion  through  each  smoke  at  visual  (0.4  -  0.7  (Jm),  near  IR  (1.060  ±  0.008  pm), 
mid-IR  (3.443  ±  0.079  pm) ,  and  far  IR  (9.750  ±  2.121  pm)  wavelengths. 
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METHOD 


Data,  in  the  form  of  time  and  CL,  and  time  and  transmittance  for  four  wave¬ 
lengths,  was  transferred  from  seven-track  tapes  to  cassette  data  tapes  for  use 
with  an  HP9830  desk  top  calculator  and  an  HP9876  plotter.  This  facilitated 
virtually  no-cost  data  analysis.  Measurements  of  CL  and  transmittance  were 
taken  at  one  to  two  second  intervals,  but  were  not  synchronized  to  each  other. 
Therefore,  the  CL  values  were  linearly  interpolated  to  the  time  of  the  recorded 
transmittance  measurement.  As  each  point  (CL  versus  transmittance)  was  plotted, 
a  normal  equation  coefficient  matrix  was  formed,  which  allowed  a  regression 
analysis  to  be  performed  on  not  only  each  trial,  but  by  adding  the  individual 
normal  matrices,  on  combined  trials  as  well. 

In  addition  to  each  entire  trial,  analysis  was  made  of  the  data  in  one  to  two 
minute  blocks.  In  this  way,  an  apparently  anomalous  series  of  points  (due  to 
saturation  below  a  certain  transmission  level,  for  instance)  could  be  avoided, 
or  analyzed  separately,  if  desired.  Ultimately,  a  subjective  determination  of 
the  best  (the  "cleanest”  typical  or  representative)  one  or  two  minute  interval 
was  selected  for  inclusion  in  the  final  analysis. 


ANALYSIS 


There  are  two  points  of  view  regarding  the  relationship  between  CL  and  trans¬ 
mittance  and,  therefore,  the  analysis  of  the  data.  Physically,  transmittance 
(T)  is  a  function  of  the  concentration  loading  (CL),  i.e.,  T  =  f(CL)  However, 
to  the  artillery  commander,  who  wants  to  know  the  CL  necessary  to  achieve  a 
certain  transmittance  level,  CL  is  considered  the  dependent  variable,  i.e., 

CL  =  f(T).  Mathematically,  of  course,  either  approach  is  valid,  but  the  least 
squares  analysis  leads  to  slightly  different  answers  in  each  case.  Results 
from  regressing  both  CL  on  T  and  T  on  CL  are  included  in  this  report  for  the 
sake  of  completeness. 

Let  an  incident  beam  of  light,  I  ,  transit  a  path  length  of  d^(m),  through  a 
medium  of  density  p(gm/m3),  with  a  mass  extinction  coefficient  k(m2/gm),  the 
sum  of  absorption  and  scattering  processes.  The  attenuation  of  the  beam  is 

dl  =  klpds, 

from  which  it  follows  that 

dI  - 

=  kpds 

or 

d  (In  I)  =  kpds. 

Integrating, 

In  I  -  In  I  *  kps, 

where  In  Iq  is  the  constant  of  integration.  Thus 

In  j-  =  kps 
o 

I_ 

is  the  fitting  equation  for  the  linear  least  squares  analysis,  with  Iq  defined 
as  the  transmittance  (in  fraction  of  incident  beam  transmitted),  and  ps(gm/m2) 
defined  as  the  concentration  loading,  also  referred  to  as  cross  wind  integrated 
concentration  (CWIC)  or  column  density.  Thus 
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In  T  =  kCL 


or 


T  = 


kCL 

e 


For  the  case  where  CL  is  considered  the  dependent  variable,  the  fitting 
equation  is 

CL  =  k'lnT. 

If  the  correlation  between  CL  and  InT  were  perfect,  k'  would  equal  r.  Since 

K. 

it  is  not,  both  k'  and  k  are  shown  in  Tables  2-4. 


In  addition  to  k'  and  k,  Tables  2-4  present  the  standard  errors  (e)  of  k'  and 
k,  and  the  standard  error  of  estimate  (Se)  of  the  dependent  variable  on  the 
independent  variable  for  each  trial  and  for  combined  ensembles  of  the  same 
smoke  type.  The  traditional  definition  of  Se  and  e  are 


Se 


r  N 

I 

i=l 

(y.  -  ax.)2 

J  l  i 

N-l 

and 


respectively,  where  y^  -  ax^  are  the  residuals  and  N  is  the  total  number  of 
points . 

For  FWP  it  may  be  inappropriate  to  compute  ensemble  averages  since  the  smoke 
types  may  be  slightly  different  for  each  FWP.  Combining  each  of  the  HC  and 
the  WP  trials,  however,  should  be  valid  because  presumably  they  are  the  same 
smoke  types. 

In  order  to  test  the  hypothesis  that  each  of  the  two  trials  for  a  given  smoke 
type  is  drawn  from  the  same  statistical  population,  a  two-sample  t  test  for  k 
and  k' ,  and  a  two-tail  F  test  on  the  e’s,  were  performed.  In  all  cases,  the 
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tests  indicate  overwhelmingly  that  the  two  trials  are  not  taken  from  the  same 
population.  This  is  interpreted  to  mean  that  either  the  smoke  types  used  are 
not  exactly  the  same  (e.g.,  trial  2-33  and  2-4  did  not  involve  the  same  HC), 
or  that  other  variables  (e.g.,  temperature,  humidity,  sun  angle,  etc.)  in  the 
experimental  process  were  not  constant  from  one  trial  to  the  other.  Additional 
results  from  subsequent  inventory  smoke  tests  may  shed  light  on  the  reason  for 
the  discrepancies  between  the  k‘s  and  k’’s.  As  it  stands,  the  potential  user 
of  this  information  is  left  with  the  option  of  choosing  a  single  trial  or  the 
ensemble  as  representative  of  the  smoke  characteristics.  Detailed  data  sur¬ 
rounding  each  trial  can  be  found  in  reference  a. 
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DISCUSSION 

Figure  1  shows  transmittance  as  a  function  of  CL  for  the  HC  and  WP  _,moke  types 
at  each  of  four  wavelengths.  For  each  smoke  and  wavelength,  two  lines  are 
drawn  representing  k  and  1 / k ’ ,  respectively.  It  can  be  seen  that  at  visual 
and  near  IR  wavelengths  HC  is  not  as  opaque  as  WP  (or  FWP) ,  but  at  longer 
wavelengths  the  reverse  is  true. 

Figures  2,  3,  4,  and  5  depict,  for  each  wavelength,  the  individual  trial's  k 
and  1/k'  (dashed  lines),  the  combined  k  and  1 / k '  (middle  solid  lines),  and  the 
standard  error  of  estimate  for  this  combined  k  and  1/k'  (outside  parallel 
solid  lines)  for  HC,  WP,  and  FWP.  Note  that  one  of  the  three  trials  of  FWP 
(DPI-005-4)  appears  anomalous  compared  to  the  other  two.  Since  this  particu¬ 
lar  trial  contributes  fewer  points  to  the  ensemble,  however,  the  combined 
result  is  weighted  towards  the  other  two  similar  trials. 

Figures  6  and  7  illustrate  the  frequency  (v)  dependence  of  combined  k  and 
1/k',  respectively.  A  relationship  with  a  correlation  coefficient  of  >0.99  is 
revealed  for  HC.  The  extinction  coefficients  are  related  to  v  for  HC  by 

k(HC)  =  -1.28  x  10"29  v2 
and 

1/k' (HC)  =  -1.44  x  10"29  v2 

The  W?  and  FWP  smoke  types  show  a  reversal  in  the  trend  of  k  and  1/k'  for  far 
IR;  also  noticeable  in  Figure  1,  the  WP  smoke  types  are  less  opaque  at  far  IR 
wavelengths  than  they  are  at  mid-IR.  Because  of  this,  no  attempt  was  made  to 
find  a  simple  relationship  between  k  (1/k')  and  v.  The  V2  dependence  of  k 
(1/k')  for  HC  is  at  this  time  unexplained  and  rather  surprising.  It  is  known 
that  Rayleigh  scattering  follows  a  v4  dependence  and  absorption  follows  a  v1 
dependence,  but  only  a  unique  combination  of  particle  size  distributions  and 
scattering  and  absorption  properties  would  result  in  the  extinction  coeffi¬ 
cient  being  proportional  to  v2.  In  fact,  the  most  surprising  aspect  of  the 
result  is  that  the  extinction  coefficient  depends  on  an  integer  power  of  the 
frequency  for  any  combination  of  size  distribution  and  extinction  properties. 
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TABLE  1.  Summary  of  Smoke  Week  I  Test  Trials 
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Number  of  submunitious  or  canisters 


TABLE  2.  Analysis  of  HC  for  CL  =  k’lnT  and  InT  =  kCL 


Trial  No.  of 

Identification  Points  k'  ±  e  Se  k  ±  e  Se 


Visual 


DPI-002-33 

67 

-0.216  ±  0.005 

0.090 

-4 . 446  ±  0.110 

0 . 409 

DPI-002-4 

75 

-0.291  ±  0.017 

0. 148 

-2.727  ±  0.161 

0.452 

Combined 

142 

-0.232  t  0.007 

0.133 

-3.826  ±  0.116 

0.539 

Near  IR 


DPI-002-33 

182 

-0.895  ±  0.012 

0.419 

-1.081  ±  0.015 

0.461 

DPI-002-4 

126 

1 

Ut 

CO 

1+ 

o 

o 

ro 

CO 

0.447 

-0.804  ±  0.020 

0.3'3 

Combined 

308 

-0.941  ±  0.013 

0.481 

-1.007  ±  0.014 

0.498 

Mid-IR 


DPI-002-33 

181 

-6.820  ±  0.154 

0.725 

-0.134  ±  0.003 

0.102 

DPI-002-4 

199 

-9.080  ±  0.179 

1.515 

-0.102  ±  0.002 

0.161 

Combined 

380 

-8.547  ±  0.133 

1.300 

-0.107  ±  0.002 

0.145 

Far  IR 


DPI-002-33 

194 

-19.116  ±  0.534 

0.874 

-0.045  i  0.001 

0.043 

DPI-002-4 

181 

-26.184  ±  0.532 

1.568 

-0.036  ±  0.001 

0.058 

Combined 

375 

-25.516  i  0.404 

1.361 

-0.037  ±  0.001 

0.053 

TABLE  3.  Analysis  of  WP  for  CL  =  k'lnT  and  InT  =  kCL 


Trial  No.  of 

Identification  Points  k'  ±  e  Se  k  ±  e 


Visual 


DPI-002-24 

116 

-0.372 

+ 

0.013 

0.249 

-2.341 

+ 

0.084 

DPI-002-34 

113 

-0.553 

+ 

0.011 

0.289 

-1.732 

+ 

0.034 

Combined 

229 

-0.493 

+ 

0.010 

0.325 

-1.854 

+ 

0.038 

Near  IR 


DPI-002-24 

180 

-1.005 

+ 

0.014 

0.361 

-0.963 

+ 

0.013 

DPI-002-34 

128 

-1.354 

+ 

0.021 

0.299 

-0.717 

+ 

0.011 

Combined 

308 

-1.082 

+ 

0.014 

0.418 

-0.880 

+ 

0.011 

Mid-IR 


DPI-002-24 

268 

-5.032 

+ 

0.057 

1.367 

-0.192 

+ 

0.001 

DPI-002-34 

224 

-5.904 

± 

0.097 

1.589 

-0.160 

+ 

0.003 

Combined 

492 

-5.311 

± 

0.054 

1.564 

-0.179 

+ 

0.002 

Far  IR 


DPI-002-24 

269 

-4.721 

+ 

0.050 

1.293 

-0 . 206 

+ 

0.002 

DPI-002-34 

221 

-5.231 

+ 

0.084 

1.552 

-0.181 

+ 

0.003 

Combined 

490 

-4.893 

+ 

0.046 

1.456 

-0.196 

+ 

0.002 

TABLE  4.  Analysis  of  FVP  for  CL  =  k'lnT  and  InT  =  kCL 


Trial  No.  of 

Identification  Points  k'  ±  e 


k  ±  e  Se 


Visual 


f  ,  V 

DPI-005-4 

(no  data  available  for 

visual  wavel 

engths) 

►  V 

DPI-005-12 

176 

-0.410  ±  0.010 

0.229 

-2.202  ±  0.055 

0.531 

!w 

J- 

DPI-005-16 

235 

-0.359  ±  0.006 

0.111 

-2.622  t  0.043 

0.299 

i 

.  < 

Combined 

411 

-0.389  ±  0.006 

0.175 

-2.257  ±  0.026 

0.430 

t 

* 

;■  t 

Near 

IR 

:  <- 
3 

DPI-005-4 

87 

-1.200  ±  0.019 

0.  182 

-0.815  ±  0.013 

0.150 

V 

DPI-005-12 

256 

-1.208  ±  0.023 

0.531 

-0.756  ±  0.015 

0.420 

DPI-005-16 

205 

-1.049  ±  0.010 

0.217 

-0.935  ±  0.009 

0,205 

:  * 

; 

Combined 

548 

-1.141  ±  0.013 

0.407 

-0.822  ±  0.009 

0 . 345 

^  1 

Mid- 

IR 

f  4 

DPI-005-4 

109 

-1.792  ±  0.062 

0.416 

-0.495  ±  0.017 

0.218 

f 

DPI-005-12 

254 

-4.421  t  0.075 

0.550 

-0.211  ±  0.004 

0.120 

DPI-005-16 

240 

-4.548  ±  0.057 

0.872 

-0.212  t  0.003 

0.188 

r- 

| 

Combined 

603 

-4.153  ±  0.053 

0.973 

-0.219  ±  0.003 

0.224 

5 

** ' 

Far 

IR 

> 

: 

DPI-005-4 

151 

-2.091  ±  0.054 

0.316 

-0.435  ±  0.011 

0.144 

DPI-005-12 

254 

-3.798  ±  0.059 

0.505 

-0.248  i  0.004 

0.129 

Y 

; 

DPI-005-16 

241 

-4.288  ±  0.047 

0.799 

-0.227  1  0.003 

0.184 

i 

Combined 

646 

-4.005  ±  0.039 

0.772 

-0.235  ±  0.002 

0.187 
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FIGURE  1.  Transmittance  vs.  Cl.  tor  comb  I  net 


FIGURE  2.  LnT  vs  CL,  Visual  Wavelengths 


FIGURE  4.  LnT  vs  CL,  Mid  Infrared 


FIGURE  6.  Extinction  Coefficient  (k)  vs  frequency  (»/) 


